Mon. Not. R. Astron. Soc. OOP. [TlfTSl C2008') Printed 24 October 2008 (MN 1^1^ style file v2.2) 



The low-mass Initial Mass Function in the young cluster 
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ABSTRACT 

NGC6611 is the massive young cluster (2— 3Myr) that ionises the Eagle Nebula. We 
present very deep photometric observations of the central region of NGC 6611 obtained 
with the Hubble Space Telescope and the following fihers: ACS/WFC F775W and 
F850LP and NIC2 FllOW and F160W, loosely equivalent to ground-based IZJH 
filters. This survey reaches down to / ~ 26 mag. We construct the Initial Mass Function 
(IMF) from ~1.5Mq well into the brown dwarf regime (down to ~O.O2M0). We 
have detected 30 — 35 brown dwarf candidates in this sample. The low-mass IMF 
is co mbined with a highe r-mass IMF constructed from the groundbased catalogue 
from lOliveira et all (|2005f ). We compare the final IMF with those of well studied star 
forming regions: we find that the IMF of NGC 6611 more closely resembles that of the 
low-mass star forming region in Taurus than that of the more massive Orion Nebula 
Cluster (ONC). We conclude that there seems to be no severe environmental effect in 
the IMF due to the proximity of the massive stars in NGC 6611. 

Key words: stars: late-type - stars: low-mass, brown dwarfs - stars: luminosity 
function, mass function - stars: pre-main-sequence - open cluster and associations: 
individual: NGC 6611. 
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1 INTRODUCTION 

The most favourable populations to study the low-mass Ini- 
tial Mass Function (IMF) are young, dynamically unevolved 
clusters, not only because they sample the entire mass spec- 
trum but also because very-low-mass stars and brown dwarfs 
are intrinsically brighter at young ages. Properties of in- 
terest in the observed IMFs of young associations are the 
higher-mass slope, the mass at which the IMF reaches a 
maximum (turn-over mass or characteristic mass) and the 
brown dwarf to star ratio. The local Universe high-mass IMF 
seems to be universal within the expected uncertainties, fol- 
lowin g a power-law with approximately the Salpe ter index 
(e.g.. iMassev. Johnson fc DeGioia-Eastwoodlll995l : for a re- 



view see 



ElmegreenI 20081 ). However, at the lower mass end 



important differences have become apparent. 

In recent years, substantial advances have been made 
in the study of the low-mass and substellar I MF, espe- 
cially in four nearby star forming regions: Taurus 
200J), IC348 jLuhman et al.ll2003l '). Chamael eon I 



Luhman 



Luhman 



2007 ) and the Orion Nebula Cluster fOlNfC. [Muench et al 



20021 : ISlesnick. Hillenbrand fc Carpenteij 12004 ). These pop- 
ulations sample different star-forming environments, from 
Taurus that represents the best studied example of star for- 
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mation at low (stellar and gas) density to the ONC that 
is a relatively extreme environment (when compared to 
these other associations), being much denser and having 
a high-mass population that extends to massive 0-stars. The 
low-mass IMF for these four regions shows important vari- 
ations, name ly in the turn-over meiss and the b rown dwarf 
to star ratio (|Luhman et al.ll2006l : lLuhmanll2007l ). This sug- 
gests a dependence of the IMF properties on the conditions 
of st ar formation. Star fo r mation theories by fragmenta- 
tion (iBate fc Bonncll"2005': 'Bonn ell et al.ll2007^ . core ejec- 
tion (e.g., R cipurth fc Clarke 2 00 ll ) and core evaporation 
(e.g., IWhitworth fc ZinneckeijbOM ) predict changes in the 
brown dwarf fraction and/or characteristic mass of the IMF 
depending on local environmental conditions (e.g., molecu- 
lar cloud density, radiation field). Constraining any environ- 
mental dependence of the IMF might help pin down which 
physical parameters are important in shaping it. 

The massive stars in the young cluster NGC 6611 
are responsible for the ionisation of the Hll region M 16, 
the Eagle Nebula. Indeed, just the 3 most massive stars 
in the cluster central area emit ~ 5 times more ionis- 
ing radiation than the central Trapezium cluster of the 
ONC. NCC6611 cluster members are distributed over a 
region of ~14arcmin radius, with a higher concentration 
in the largely uno bscured 4arcmin radius central area 
l|Belikov et al] I2OO0I ). Recent distance determinations, de- 
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rived using spectroscopic p arallaxes JPufton at al 



and main-sequence t urnoff (|Bonatto. Santos fc Bical 



2006) 



200e 



iGuarcello etaP l2007[ l. favour values around 1.8 kpc. The 
cluster contains a large number of massive stars as well 
as a large population of pre-main sequenc e (PMS) stars 
ijHillenbrand et al.lll993l : l01iveira et al.ll2005l ). The total ob- 
served mass of NGC6611 is at least a factor ^ 2 large r 
than that of the ONC (|Bonatto. Santos fc Bical |2006| ). 
This cluster is very young, with a i i age of 2—3 Myr 
iHillenbrand et al.lll993i : iBelikov et allbOOOl ). A consider- 
able age spread (< 1 —6 Myr) has been rep orted fo r sources 
in the Eagle Nebula (iHillenbrand et al.l[l99 3: Bcli kov et all 
I2OO0I : iDufton et al.ll2006l) but llndebetouw et all (|2007l ) find 
no clear evidence of age gradients throughout the region. 
A rich low-mass (M > 0.5 M ) PMS population was identi- 
fied bv lOliveira enn(|2005l ) concentrated towards the most 
massive NGC6611 cluster members. Many of these young 
stars retain dusty circumstellar discs, identified by L-band 
excesses. For a co mprehen s ive re view on NGC6611 and the 
Eagle Nebula see lOliveiral (|2008l ). 

In this paper we discuss new Hubble Space Telescope 
(HST) observations of the central area of NGC 6611 at opti- 
cal and near-IR wavelengths. These new observations are the 
deepest to date in this region, reaching well into the brown 
dwarf regime. We construct the low-mass IMF from inter- 
mediate to substellar masses by combinin g the HST dataset 
with ground-based photometry from Oliveira et al.l (|2005h . 
This manuscript is organised as follows. Firstly we describe 
the HST data and its reduction process and the photometry 
and its calibration. This is followed by the PMS identifica- 
tion, reddening determinations and completeness and con- 
tamination corrections, both for the HST data and ground- 
based data. Finally, we construct the IMF for NGC 6611 and 
compare its main features (slope, characteristic mass and 
brown dwarf to star ratio) to other well studied young clus- 
ters and associations. We discuss the implications in terms 
of the environmental impact on the IMF and star and brown 
dwarf formation scenarios. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 HST observations 

The Hubble Space Telescope images of the core of NGC 6611 
described here were obtained as part of the cycle 14 pro- 
gramme number 10533 (P.I. J.M. Oliveira). Optical images 
were obtained with the Wide Field C hannel (WFC) of the 
Advanced Camera for Surveys (ACS, iBoffi et al.ll2007D . A 
single ACS/WFC field was observed through the F775W 
and F850LP filters (roughly corresponding to ground-based 
I and Z). With a ~ 0.05"pixel scale and two 2K x 4K detec- 
tors, the WFC field-of-view is ~ 3.4' x 3.4'. The system gain 
was set to 2 to completely sample the full-well-depth (as 
opposed to ~ 75% for gain 1), resulting in only a modest in- 
crease in readnoise and a significant increase in the dynamic 
range of the observations. The total exposure time per fil- 
ter is 2 000 s, split into 4 dithered exposures of 500 s. This 
corresponds to one HST orbit per filter. The observations 
were designed to probe well into the brown dwarf regime, 
reaching down for the first time in such a distant cluster to 
~O.O2M0. 



Figure 1. Top: true colour image of NGC 6611 showing the 
ACS/WFC field (I-band in blue and Z-band in red). The field 
centre is at 0=18:18:41.72 & (5=-13:47:35.0. The Eagle Nebula 
pillars are located to the south-east of this field, approximately 
3' away. Bottom: true colour image showing the NICMOS mosaic 
(J-band in blue and H-band in red). 



Near-Infrared images were obtained with the Near In- 
frared Camera and Multi-Object Spectrometer (NICMOS, 
iBarket et all |2007| ) onboard HST. The camera used was 
NIC2 with a 19.2" x 19.2" field-of-view and 0.075" pixels. 
We observed 25 NIC2 fields to form a mosaic covering 
^ 2.5 arcmin'^. The NIC2 field positions were chosen in such 
a way as to avoid the bright stars in the field (Fig.[T} . The fil- 
ters used are FllOW and F160W. These correspond loosely 
to ground-based J and H filters, but the reader should note 
that the FllOW is much broader than any other commonly 
used J-band filters (0.8— 1.4 /xm). The total exposure time 
per field per filter is 128 s, split into 4 dithered exposures. To 
observe all 25 fields in both filters 5 HST orbits were needed. 
The area covered by the HST observations in NGC 6611 is 
shown in Fig.[T] 
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2.2 ACS/WFC calibrations and PSF photometry 

The ACS reduction pipeline consists of two main stages: the 
first corrects for instrumental effects and produces calibrated 
data products including appropriate photometry keywords 
(CALACS); the second corrects for the significant geometric 
distortion in WFC images, performs cosmic ray rejection and 
combines asso ciated dithered data u sing drizzle techniques 
(MultiDrizzle, IPavlovskv et al.l 120061 ). For our programme, 
the HST pipeline results were a single astrometrically and 
photometrically calibrated image for each filter ("drz im- 
ages" in units of e~s~^). The photometric zeropoints are 
25.256 and 24.326 mag r espectively for I a nd Z (VEGA- 
MAG magnitude system, ISirianni et al. I l2005h . Point Spread 
Function (PSF) photometry was performed on these pipeline 
images as described below. 

PSF fitting packages use an optimal weighting scheme 
that depends on the true counts per pixel as well as readout 
noise and gain. To properly compute the charge transfer 
efficiency (CTE) correction (see below) it is also necessary 
to have the total number of electrons instead of electron rate. 
Therefore, before applying the photometric algorithms, the 
drz images were multiplied by the total exposure time and 
the average background sky level was added back to the 
images (see extensive discussion bv lSirianni et al.ll2005l ). 

Stell ar photometry was obtained us ing the DAO PHOT 
package l|Stetsonlll98"% within P^j;Ajfl (|Davislll994 ), mod- 
elling the spatially variable PSF independently for each fil- 
ter. Standard object detection and aperture photometry al- 
gorithms were used to detect and compute preliminary mag- 
nitudes. PSFs were constructed using a large number of stars 
well distributed across the images (respectively 134 and 160 
stars for the I- and Z-bands). Stars were chosen to be bright 
so that the PSF in the core and wings could be well de- 
scribed but not brighter than ~18.5mag; above this limit 
there are clear non-linearity effects and the PSF profiles 
of those stars no longer match the profiles of fainter stars. 
An elliptical Moffat function was chosen for the functional 
form of the analytical component of the PSF model and this 
model was allowed to vary quadratically with position in 
the image. Once appropriate PSF models were computed, 
the ALLSTAR function was used to compute photometry 
for all the stars in the frame; approximately 16 000 objects 
have measured PSF photometry in each frame. As can be 
seen in Fig.[T] the ACS images contain a number of bright 
stars that saturate the exposures and cause leakage artifacts. 
These artifacts cause a large number of spurious detections 
distributed mostly in halos around the bright stars and on 
spikes along affected columns and/or rows. 

The PSF photometry output files contain two output 
parameters that together with the magnitude error can be 
used to filter out spurious detections. The sharpness param- 
eter estimates the intrinsic angular size of the object: stellar 
objects should have a value around zero (with a larger scat- 
ter for fainter magnitudes) , while large negative and positive 
values identify respectively cosmic rays and other blemishes 
and resolved galaxies and blended objects. The param- 
eter measures the goodness-of-the-fit and should be around 



unity with no noticeable trend with magnitude. As men- 
tioned above, bright but unsaturated stars are affected by 
non-linearity effects so the measured is somewhat larger 
for these objects. A very good description on how imposing 
cuts to these parameters is used to eliminate spurious de - 
tections in catalogues can be found in lReikuba et al.l (|2005h . 
The photometric errors, sharpness and x^ parameters were 
thus used to refine each photometry catalogue and then the / 
and Z detections were cross-correlated. As both filters were 
observed in the same visit the positional drifts between the 
images are negligible, so the catalogues were matched in X 
and Y pixel positions before astrometry was performed: 3462 
objects were detected in the combined IZ catalogue. 

The ACS CCD detectors suffer from scattered light at 
long wavelengths (from about 7 500 A) and the fraction of 
integrated light in the scattered light halo increases as a 
function of wavelength. The amount of light contained in 
the halo will also be larger the redder the stellar source, 
i.e. the shape of the observed PSF depends on the colour 
of the star, something that PSF fitting algorithms cannot 
account for. The best way to mitigate for this is to compute 
colour-dependent aperture corrections. Therefore, we were 
only able to calibrate in this way objects that are detected 
in both the I- and Z-band filters. The halo effect and the 
prescription on how to correct for it are described in detail in 
ISirianni et al.l (|2005h and we provide here only a summary. 

A useful quantity to help estimate colour-dependent 
aperture corrections is the effective wavelength Aeff that rep- 
resents the mean wavelength of the detected photons. This 
source-dependent passband parameter can be used to quan- 
tify changes in the stellar light encircled in a radius r as a 
function of the stellar colo ur. Once Aeg i s det ermined for 
each source and each filter, ISirianni et al.l (|2005l ) provides a 
tabulated relation between Aoff and aperture correction for 
several aperture radii. 

To define a relation between AoH and observable in- 
stru mental colour we u sed the function calcphot in syn- 
phot (|Laidler et al.ll2005l ) and a large number of stellar tem- 
plate spectra. We used the following spectral atlases: the 
Bruzual-Persson-Gunn-Stryker atlas available within syn- 
phot (spectral types to early M) and the atlas of M, L and 
T type dwarfs com piled by Sandy Leggetl[f| with spectra 
from several auth ors l|Leggett et al.|[200Gl : lKnapp et al.ll2004l : 
IChiu et al.ll2006l ). In total over 300 spectra were used from 
early-type stars to T-dwarfs. With calcphot, we computed in- 
strumental I—Z colours, Aeff, I and XeB,z to co nstruct a linear 
relation equivalent to that shown in Fig. 11 of lSirianni et al.l 
(l2005h . For the objects in our combined IZ catalogue, we 
measured their instrumental colour, interpolated to esti- 
mate Aoff, I and Aoff,z, and interpolated Sirianni's Table 6 
to obtain colour-dependent aperture corrections. Computed 
aperture corrections are in the ranges 0.22—0.25 mag and 
0.42—0.55 mag, respectively for / and Z. 

Another effect that needs to be corrected for is the 
position- and time-depen dent decrease in the CCD's charge 
transfer efliciency (CTE. lRiess fc Mackir2004l '). A photomet- 
ric correction can be applied that depends on the stellar flux 
and sky levels, on the Y position in the chip and on the time 
interval since the start of ACS operations. Using the pre- 



^ PyRAF is a command language for running IRAF tasks that 
is based on the Python scripting language. 



^ http:/ /www. jach.hawaii.edu/~skl/LTdata. html 
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scription and coefficients described in lRiess fc Maclj (|2004l ) 
we computed magnitude corrections of up to 0.05 mag and 
0.08 mag respectively for I and Z. 

Once all these photometric corrections were applied, as- 
trometry was performed using the IRAF task xy2rd. As- 
trometry was performed individually in the / and Z images 
and the positions were averaged. No absolute astrometry 
was performed at this stage. The final IZ catalogue reaches 
down to J = 26 mag and Z = 2A mag, with photometric 
errors '--^ 0.12 mag and ~ 0.08 mag respectively. 



2.3 NIC2 mosaics, calibrations and aperture 
photometry 

The NICMOS pipeline consists of two main stages: the first 
corrects for instrumental effects (CALNICA) and the sec- 
ond combines images obtained in a dither pattern and sub- 
tracts the sky background from the images (CALNICB). 
Due to the way in which the observations were designed, the 
pipeline combines the 4 dithered exposures for each field but 
it does not combine the 25 fields in a single mosaic. Crucially 
the pipeline does not correct for the geometric distortion, 
the largest component of that distortion being a large dif- 
ference in the X and Y platescales. The best way to correct 
for these distortions and simultaneously construct complete 
mosaics for each of the J and H filters is to use the rou- 
tine multidrizzle in PyRAF on the output frames from the 
CALNICA step, i.e on the 25 x 4 images per filter. This rou- 
tine works in a very similar way to the second pipeline step 
for the ACS images. Mosaics were constructed for J and H 
with a final pixel scale of 0.05", that conveniently matches 
the WFC pixel scale. 

Before the images were mosaiced they were corrected for 
the NIC MOS count-rat e- and wavelength-dependent non- 
linearity (jde^ Jongj 120061 ). This correction is not applied by 
the pipeline reduction process but a separate Python rou- 
tine is availablfl The non-linearity correction amounts to 
0.063 mag offset per dex change in incident flux for J and 
0.029 mag per dex in H. The routine changes the count 
rates in the output images and therefore new zero-points 
are needed: 24.561 and 23.863 mag respectively for J and H 
(VEGAMAG magnitude) for images in units of e~s~^. 

We initially attempted to perform PSF photometry on 
the J and H images as well. However, NIC2 undersam- 
ples the PSF a t these wavelengt hs and the PSF profile is 
rather complex l|Krist et al.lll998l ): as a result we found that 
the quality of the PSF fits was not good, particularly for 
the H-band. As our images are not crowded, we opted to 
use aperture photometry instead. Approximately 5 000 and 
10 000 sources were detected in J and H respectively; a large 
number of these are spurious detections (e.g., artifacts near 
bright stars and at the edge of individual frames etc) that 
can be removed by combining these catalogues with the IZ 
catalogue. Aperture corrections were computed using 60—80 
unsaturated bright stars. Astrometry was performed on the 
J and H catalogues as described above for the ACS im- 
ages. Possible deficiencies in using aperture photometry are 



http:/ /www.stsci.edu/hst/nicmos/performance/anomalies/ 
nonlinearity. html 



corrected for with a careful computation of completeness 
corrections described in Section 3.1.3. 

2.4 Ground-based complementary data 

To complement our HST catalogue at the brighter end (the 
HST images saturate for PMS objects of just over a solar 
mass), we used the optical and IR catalogue described in 
lOliveira et all Hooi), that covers a much larger area in the 
central region of NGC6611. For a complete description of 
the observations and data reduction the reader should refer 
to that article. Their catalogue includes I and Z Cousins 
magnitudes and JHK photometr y in the MKO-IR system 
(|Tokunaga. Simons fc Vaccall2002l ). 

As noted above the HST filters are very different from 
groundbased fllter sets and it is not trivial to convert be- 
tween those systems as significant colour terms are expected 
(especially for cool PMS objects) due to the very different 
filter bandwidths. Even though there are a number of objects 
observed in both HST and groundbased filter sets, these are 
few and crucially they are at the bright/blue end of the HST 
data and therefore the colour range is not sufficient to de- 
rive reliable magnitude conversions. Furthermore we have 
to use two different sets of PMS isochrones, as no single 
set covers the whole mass range we are addressing in this 
manuscript, from several solar masses down to the substel- 
lar regime. Therefore we have made no attempt to merge 
these two photometric catalogues. 

The groundbased catalogue reaches down to ~0.5Mq 
and thus provides a good overlap in PMS stellar masses 
with the HST catalogue. As described below we construct 
cluster IMFs over different mass ranges (over the same area 
and using methods that are similar whenever possible) and 
combine them for the total IMF analysis. 

3 COLOUR-MAGNITUDE AND 

COLOUR-COLOUR DIAGRAM ANALYSIS 

3.1 Low-mass stars and brown dwarfs 

The IZ, J and H HST catalogues were combined based on 
their astrometry. As can be seen in Fig. [1] the J and H im- 
ages cover roughly a quarter of the I and Z images. In the 
overlapping region 729 objects have IZJH photometry — / 
band magnitudes are the limiting ones. This is the HST cata- 
logue from which pre- main-sequence stars are going to be se- 
lected based on several colour-magnitude diagrams (Fig. [2]). 

3.1.1 Identification of PMS stars and membership 
selection 

Fig. [2] shows the colour-magiutude diagrams used to iden- 
tify the pre-main-sequence (PMS) population in NGC6611. 
These young stars are cooler than field stars of the same ap- 
parent magnitude and thus can be easily identified in such 
diagrams. Selection criteria are shown in the figure, based on 
the observed I magnitude and I — Z, I — J, and J — H colours 
and the position of theoretical isochrones (|Baraffe et al.l 
Il998l l. Using these criteria 290 stars are identified as PMS 
candidates and as such, candidate low-mass members of 
NGC6611 (Table[T] complete in electronic form only). A 
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Figure 2. Colour-magnitude diagrams of the HST survey area. The straight hues indicate the regions for selection of PMS candidates: 
open circles are all objects with IZJH photometry while filled circles mark the PMS candidates. Note the large number of heavily 
obscured stars towards the Galactic centre with J — H ~ 2 — 3 mag. PMS isochrones for 2 Myr (full line) and 3 Myr (dashed line) from 
19981), assuming average extinction (Sect. 3.1.2) and a distance of 1.8 kpc are also shown. Red lines are isochrones from 
1 2008h for distance and extinction consist ent with an evolved population in the Scutum-Crux spiral arm (see text). Blue 



Baraffe et al. 



Marigo et al 



squares and stars are Upper Scorpius llKraus et al .|2005^ and Taurus l lKraus et al. l200d) very low-mass and brown dwarf cluster members, 
of ages respectively ~5Myr and 1—2 Myr, observed with the same filters and placed at the same distance and average extinction as 
NGC 6611. 



Table 1. NGC 6611 cluster candidates selected from the HST photometric catalogue. Absolute astrometr y was performed on this 
catalogue using the positions of ~ 40 bright stars that are also detected in the catalogue of lOliveira et al.l l l2005^ . that was astrometrically 
calibrated against 2MASS. The precision in the positions is better than 0.04" in each coordinate. Columns provide the magnitudes and 
errors for /, Z, J and H. E(B — y)2Myr and E{B — V)3Myr a^re the extinction determinations using the 2 or 3 Myr isochrones. Based 
on these determinations, candidates are classified as members (0.45 < E{B — V) < l.Omag) or non-members (see Section 3.1.2). M2Myr 
and MsMyr are the masses of PMS objects, assuming either an age of 2 or 3 Myr. The complete table is available in electronic form. 

ra dec IZJH E(B — V)2Myr E(B — V)3iy[y,- status M2Myr A'/sMyr 

") °) (mag) (mag) (mag) (mag) (mag) (mag) (Mq) (Mq) 

18:18:46.09 -13:43:01.3 17.842 ±0.021 16.754 ±0.011 15.674 ±0.002 14.354 ±0.002 0.54 ±0.08 0.55 ±0.06 member 0.678 0.750 

18:18:46.06 -13:43:02.3 19.032 ±0.012 17.893 ±0.015 16.846 ±0.004 15.500 ±0.003 0.58 ±0.06 0.58 ±0.05 member 0.402 0.457 

18:18:46.02 -13:43:05.8 20.934 ±0.014 19.334 ±0.014 18.009 ±0.008 16.708 ±0.006 0.63 ±0.05 0.64 ±0.05 member 0.132 0.164 

18:18:45.75 -13:42:47.1 20.032 ±0.026 18.622 ±0.015 17.468 ±0.005 16.047 ±0.004 0.71 ±0.05 0.78 ±0.05 member 0.270 0.334 

18:18:45.71 -13:42:44.1 19.473 ±0.019 18.192 ±0.017 17.099 ±0.004 15.975 ±0.004 0.31 ±0.05 0.34 ±0.05 non member 

18:18:45.71 -13:42:47.2 19.717 ±0.019 18.354 ±0.016 17.223 ±0.004 15.861 ±0.004 0.62 ±0.05 0.68 ±0.05 member 0.284 0.350 



census of circumstellar discs e xists for the brighter 
part of this sample (J < 1 9 mag, Oliveira et al I I2OO5I : 
see also llndebetouw et al. 2007 ^ but we do not take 
this information into account in order not to bias 
our sample selection in any way. Furthermore, no 
disc survey exists for the larger, fainter part of our 
sample. 

We also show in the figure two populations of spec- 
troscopically confirmed young very low-mass stars and 
brown-dwarfs, with ages that nicely bracket the age 



of the PMS population in NGC 6611 and that have 
been observed with the same filterfl Taurus (1-2 Myr; 
iKraus. White fc Hillenbrand! |2006|) and U pper Scorpius 
(~5Mvr: lKraus. White fc Hillenbrandll2005h . The position 



* These objects were observed with the same filters but with the 
HST ACS/HRC camera; very small colour-dependent transfor- 
mations are applied to convert the ACS/HRC measurements to 
magnitudes as would be observed with the ACS/WFC camera. 
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of these objects in the diagram shows that our selection re- 
gion is generous enough to include a PMS population of age 
< 5 Myr; our selection therefore should not exclude any PMS 
objects. 

The foreground loci in the colour-magnitude diagrams 
are sparsely populated so contamination by foreground 
dwarfs is not large; this was expected since the area we are 
analysing is relatively small. The cuts shown in the figure, 
in particular the cut m J — H colour, are also very effec- 
tive in separating the PMS population from the large num- 
ber of heavily reddened stars towards the Galactic centre. 
These stars delineate a well-defined sequence with I — Z, 
I — J and J — H colours in the narrow ranges 2—3, 4.5—6 
and 2— 2. 5 mag respectively. Our HST survey is so deep that 
in order to understand the nature of this heavily reddened 
population we need to consider what lies behind the cluster 
towards the Galactic centre. 

The line-of-sight towards the Eagle Nebula encounters 
the S cutum-Crux spiral arm at a distance of < 4kpc (|Valleel 
I2OO8I ). Recent three-dimensional Gal actic insterstellar ex - 
tinction maps (based on 2MASS data. [Marshall et aLlbooel ) 
show that behind NGC6611 there is a marked extinction 
jump consistent with the location of this spiral arm: between 
2.8—3.5 kpc from the Sun, extinction as measured by Ak in- 
creases from ~ 0.6 to ~ 1.4 mag. This suggests that the dis- 
tinct clump of reddened stars could be part of the Scutum- 
Crux spiral arm struct ure. To test this hyp othesis we use 
evolutionary tracks from lMarigo et al] l|2008l ) for a 5 Gyr-old 
population and solar metallicity, computed for the HST ACS 
and NICMOS filter sets. Two such tracks are shown in Fig. [21 
the bluest for a distance of 2 500 kpc and Av — 13 mag and 
the reddest for 3 500 kpc and Ay = 17 mag. These tracks en- 
compass the observed population very well. Therefore, even 
though the distance and extinction estimates are very ap- 
proximate, the population with redder colours is entirely 
consistent with what we would expect for a population in 
that spiral arm, at a range of distances and extinctions con- 
sistent with the 2MASS extinction profile. This population, 
that would include both main-sequence and evolved stars, 
seems to be an important source of contamination. However, 
Fig. [5] shows that these stars are so much redder than gen- 
uine cluster members that they are effectively excluded by 
the J ~ H cut we apply for PMS selection. Other sources of 
contamination are discussed further in Section 3.1.4. 



3.1.2 Cluster reddening and de-reddening of PMS stars 

The first step to analyse a PMS population in a region like 
NGC6611 is to determine intrinsic magnitudes, i.e to cor- 
rect for reddening. Many authors have stu died the extinction 
properties towards the Eagle N e bula (e.g., Hillenbrand et al] 
ll993l : lBelikov et al.lll999l . i2000l : iDufton et al.ll200^ : the op- 
tical extinction (^v) is variable and the values of Ry (the 
ratio of total to selective extinction), estimated on a star by 
star basis, are larger than the normal interstellar medium 
value. Reported Ry values are in the ra nge 3.5—4.8 (typ- 
ical value ~ 3.75, iHillenbrand et ah ||1993^ while E{B - V) 
values vary between 0.5— 1.1 mag. Towards the NGC6611 
cluster extinction is however at the lower en d of this range 
(IBelikov et al. l ll999l : llndebetouw et al.ll2007l ). 

The way in which we designed the NICMOS obser- 
vations (i.e. avoiding the bright stars) means that we do 




right ascension 

Figure 3. Schematic diagram showing our survey area and pub- 
lished measurem ents of E(B — V) t owar ds O and B s t ars in 
NGC6611, from iDufton et al.l 1I2OO6I) and iBelikov et al.l lll999l . 
l2000t) . The typical value for the survey area is E{B — V) ^ 
0.70 mag, with measurements in the range 0.45 to 1 mag. 



not have independent extinction determinations towards 
the survey area. However, as can be seen from Fig.[3l that 
area is surrounded by many stars with reddening measure- 
ments. Thus a typical E{B — V) towards the survey area 
is of the order of 0.70 mag with an observed spread of 
~ 0.45—1.0 mag, corresponding to Ay ~ 2.6 mag and a range 
of ~ 1.7— 3. 75 mag — we adopt Ry = 3.75 for cluster mem- 
bers. Thus, extinction is relatively low in the central cluster 
area but it is also variable even at small spatial scales. 

Accordingly we have opted to compute extinction to- 
wards individual PMS candidates using the two colour- 
colour diagrams sh own in Fig. 2] Also shown are the 
iBaraffe et al.l (1 19981 ) PMS isochrones (these authors have 
produced absolute magnitudes for the 4 HST filters used in 
this analysis) . Using the stellar spectral atlas and procedure 
described previously, we compute HST colours for a large 
number of dwarfs in order to construct a dwarf locus down 
to spectral type M8— M9. Reddening vectors are computed 
us ing the same template spec t ra an d the reddening curve 
of ICardelh. Clavton fc MathisI (|l989l ') adopting Ry =3.75: 
the spectra are reddened with a range of E{B — V) val- 
ues and synphot is used to estimate the extincted colours 
through the HST filters. The resulting reddening slopes are 
A{I) = 2.435 X E{B - V), A{Z) = 1.86 x E{B - V), 
A(J) = 1.30 X E{B - V) and A{H) = 0.76 x E{B - V). 

PMS stars have lower surface gravities than main- 
sequence dwarf stars, and thus populate colour-colour di- 
agrams in the region between dwa rfs and giants (see also 
lOliveira. Jeffries fc van Loon^ 120041 ). As can be seen from 
Fig.|4] there are sizable differences between the isochrones 
and the determined dwarf loci. We did not have many late-M 
dwarf spectra available to constrain the reddest loci colours, 
and we do not see evidence in our dataset of the sharp turn 
to redder ,J — H colours that seems to occur for these spec- 
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Figure 4. Colour-colour dia grams for the PMS candidates selected from the HST photometric catalogue. The full and dashed blue 
lines are lBaraffe et ahl lll998l) isochrones, respectively for 2 and 3Myr. The dotted lines represented the reddening band derived for the 
HST filters. Dwarfs from an archival spectral atlas (see text) were used to empirically derive the dwarf loci (red lines, to spectral type 
M8— M9). Filled circles are the PMS objects with final reddening determination between 0.45 and 1.0 mag (see Fig.[5j. 



tral types: the variation of the J — H colour for M-dwarfs is 
~ 0.4 mag, while for our PMS candidates the observed J — H 
spread is ^0.3 mag, even with possible variable extinction 
effects. We have therefore decided to de-redden the PMS 
candidates using the set of PMS isochrones that are being 
used for the mass determination (|Baraffe et al.|[l998l ). 

Fig. |4] shows that the Baraffe isochrones in the colour- 
colour diagrams are somewhat age-dependent. The depen- 
dence is small and only affects stars in a small range of ob- 
served colours. Nevertheless, we adopt two values for the age 
of NGC 6611 , 2 and 3Myr, consisten t with the higher mas s 
population l|HiUenbrand et al.l Il993l : iBelikov et all l200d ). 
We will keep these calculations completely independent, i.e 
we obtain independent IMFs for these two assumed ages in 
order to isolate the effect of fixing the age of the population. 

Using the reddening vector and the isochrones, E{B — 
V) is computed for each star in each diagram; when avail- 
able the two estimated E{B — V) values are averaged. Error 
bars are computed for these estimates, including contribu- 
tions from photometric errors in the observed colours and 
the uncertainty introduced by the choice of the reddening 
vector — the reddening vector slope changes somewhat de- 
pending on the spectral energy distribution of the star, see 
previous discussion on AcS. 

A feature in these diagrams are the 32 stars that can- 
not be de-reddened because they do not have a reddening 
vector intersection with the isochrones. Three of these ob- 
jects are too red for the available theoretical colours. The 
remaining 29 objects have observed J — H colours that seem 
too red for their observed I — J and I — Z colours and 
thus cannot be de-reddened onto the PMS loci. Closer in- 
spection of Fig.|4] suggests that the majority of those stars 
would be too reddened to be considered bona-fide clus- 
ter members. These stars are relatively bright thus if they 
were NGC 6611 members they should have an intrinsic PMS 



(J — _ff)o < 0.9 mag. If we de-redden these objects using 
only their J — H colour to a conservative intrinsic colour 
of (J — H)o =0.9 mag, the estimated E{B — V) for most 
objects would be in the range 1.0— 1.7 mag with only three 
objects having E{B — V)< 1.0 mag. Therefore of these 29 
objects all but 3 would be rejected for having reddening too 
large when compared to the reddening distribution of more 
reliable cluster members (see Fig.[S] and discussion on red- 
dening cuts below). It is likely that these objects are actually 
behind the cluster and not PMS objects. 

Some earlier type stars fall in the region of the diagrams 
where the de-reddening procedure has two solutions. For 
these stars we computed two extinction solutions, i.e. one to 
the "late-type branch" (with 0.9 < J — H < 1.0 mag) and one 
to the "early-type branch" (with 0.6 < J—H < 0.9 mag). We 
find that the E{B — V) values determined to the "late-type 
branch" are consistent with the E{B — V) cluster distribu- 
tion (Fig. [5} . We note that the extinction determinations for 
these objects have larger error bars reflecting the uncertain- 
ties of the procedure. 

Fig.[S] shows the histograms of E{B — V) values de- 
termined towards PMS candidates. The mean of the dis- 
tribution is consistent with the average E{B — V) mea- 
sured towards massive cluster members (Fig.O. Its stan- 
dard deviation is of the order of 0.15 mag; when compared 
to the typical E{B — V) error of 0.07 mag this again sug- 
gests that local variations of extinction are real, as can be 
seen from the images themselves and Fig. [3] The range of 
E{B - V) values in Fig.[5] (0 < E(B - V) < 1.5 mag) 
is larger than the range measured towards massive cluster 
members (0.45 < E{B — V) < 1.0 mag), suggesting that 
the PMS candidate sample includes some foreground and 
background contaminants. We use these histograms to re- 
fine cluster membership: objects that have E{B — V) values 
in the range 0.45— 1.0 mag (the same as the massive clus- 
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Figure 5. Histograms of the E{B — V) determinations, using 
the 2 Myr (full line) and 3 Myr (dotted Une) isochrones. These 
determinations are consistent with the average E(B — V) towards 
massive cluster members (E{B — V) =0.7 mag). Vertical dashed 
lines show the cuts used to refine membership (see text). 



ter members) are kept in the PMS sample. Our survey is 
so deep that reddened stars behind the cluster can contami- 
nate the PMS selection (Section 3.1.1). The upper reddening 
cut is essential to exclude as many of those objects as pos- 
sible (see Section 3.1.4); the lower cut eliminates some fore- 
ground dwarfs. We note that, so long as reddening is not 
mass dependent and the intrinsic loci and reddening vectors 
are not in error in a mass-dependent way, then reddening 
cuts should not affect the derived shape of the IMF, even if 
t he y were too stringe nt . 

lOliveira et alj (j2005) determined that approxi- 
mately 55% of the objects brighter than J = 19 mag 
show an excess in the L'-band indicative of the 
presence of a circumstellar disc, while only ~ 18% 
of those objects show such an excess in the K- 
band. IZJH magnitudes are generally much less 
sensitive to disc excesses (e.g.. ' Oliveira et al.|[2005l ; 
[Robitaille et al. 2007) and therefore this should not 
significantly affect our reddening determinations. 
Nevertheless, our reddening cuts (see above) are 
generous enough to accommodate uncertainties in 
the reddening determination. 

The final PMS sample includes 208 objects. Once ex- 
tinction is computed, intrinsic 1-band magnitudes are calcu- 
lated, taking care to propagate the errors introduced by the 
extinction calculation. Then assuming either an age of 2 or 
3 Myr and adopting a distance of 1.8 kpc, masses are com- 
puted f or each object using t heir intrinsic 1-band magnitudes 
and the lBarafFe et al] (| 19981 ) isochrones. For this sample, es- 
timated masses range from 0.02 to 1.5 Mq (Table[T)). In the 
next sections we deal with two crucial issues in IMF analysis: 
completeness and contamination. 



3.1.3 Completeness correction 

Many bright stars are saturated in the HST images and cre- 
ate image artifacts that can prevent the detection of fainter 
stars nearby. We also need to characterise the reduced de- 
tection sensitivity at the fainter end. In order to address the 
issue of completeness of the photometric catalogues we con- 
struct fake star simulations: replicas of PMS stars (i.e. stars 
with true PMS magnitudes in I, Z, J and H), randomly 
selected from the observed 1-band luminosity function, are 
introduced into the 4 HST images, at random positions. For 
simplicity in generating the randomised star positions, this 
experiment was performed in the central 1.2 arcmin'^ area of 
the survey. One hundred sets of fake stars were created, each 
with 100 randomly selected PMS replicas. In each simulation 
the number of objects in the central region was increased by 
about 1/4, thus not changing the crowding properties of the 
field in any significant way. 

In order to introduce the replica PMS stars in the im- 
ages we use the PSFs determined for each filter and the 
addstar procedure in IRAF. As described in Section 2.3, for 
the J and H band we computed PSFs for the mosaics but 
ultimately opted to use aperture photometry in the final cat- 
alogue. We tested the procedure of adding stars to the mo- 
saics using the computed PSFs and then performing aper- 
ture photometry on those stars and it works well: stars are 
recovered at a similar rate as for the 1 and Z images and 
their magnitudes are preserved. 

The complete data reduction procedure, exactly as de- 
scribed in the previous sections, was repeated for each of 
the 100 simulations, including (PSF or aperture) photome- 
try, aperture and linearity corrections, astrometry, catalogue 
matching etc and PMS selection. Finally the input and re- 
covered 1-band luminosity functions are compared: of the 
10 000 stars simulated 9 301 were recovered. The recovered 
stars are not only detected in /, Z, J and H but they also 
have colours consistent with PMS selection, i.e. of the PMS 
stars simulated approximately 93% are recovered and se- 
lected as PMS objects. This is what we mean by "recovery 
rate" . 

Fig-E] shows the recovery rate as a function of 1-band 
magnitude, as well as the histograms for the simulated and 
recovered 1-band magnitudes. There is a slight magni- 
tude creep, i.e. more stars are recovered in some 
of the brig hter magnitude bin s than were originally 
simulated (^ Navlor et al.l 12002*). due either to statis- 
tical uncertainty or to th e presence of another star 
nearby (see also IStetsonI [l99ll ). As at the brighter 
end no stars are lost, this can result in a recovery 
rate larger than 1 for some bins. At the fainter end 
the recovery rate is as low as 35%. We should point 
out that in the combined IZJH catalogue completeness ef- 
fects are mostly introduced by the I and Z catalogues: the 
J- and H-band images are deeper and are less affected by 
saturation artifacts. To construct the completeness function 
we fitted a second-degree polynomial to the recovery rate 
for I < 18 mag; for brighter magnitudes the recovery is con- 
stant at ~ 1. For each object in the PMS sample a weight is 
assigned equal to the inverse of the completeness correction 
at its observed 1-band magnitude; these weights are used to 
compute completeness-corrected IMFs. 
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Figure 6. Completeness correction versus I-band magnitude. The 
liistogram (with poissonian error bars) shows the recovery rate as 
a function of magnitude (see text). The full line is a second-degree 
polynomial used to represent the completeness correction as a 
function of I-band magnitude. The inset shows the histograms of 
simulated and recovered I-band magnitudes. 



3.1.^ Field star contamination 

To make a quantitative estimate of the efltect of contami- 
nation in our catalogues we follow the procedure first out- 
lined by Jeffries et al. (2004). Firstly, a field star popula- 
tion is simulated using an I-band luminosity function for 
dwarfs in the Galactic disc. We ado pt the disc M dwar f 
I-band luminosity function $(M/) of IZheng et al.1 (|2004 ). 
complemente d at the bri ghter end with <l?(Mv) compiled 
by [Reid fc HawlevI ll2000l) . where $(Mv) is converted to 
$(M/) using averag e V — I dwarf colours l|Leggettl 1 19921 : 
iLeggett et al.l 1 19961 ). Using this luminosity function (LF 
in units of stars pc~^ mag~^) as a weighting function, ab- 
solute I-magnitudes and intrinsic colours were randomly 
assigned to a very large number of simulated stars. The 
spatial density of stars is described by A'oexp( — (i? — 
Rq) /ha) eyip{—Z /hz), where Z = \Zq -|- dsin(6)| is the 
height above the Galactic plane, Zq — -1-27 pc is the height 
of the Sun above the plane, d is the heliocentric distance, b 
is the Galactic latitude, hz ~ 270 pc is the scale height, R 
is galactocentric distance, hn = 2 250 pc is the scale length 
and A^o is the space density in the plane ( Z = 0) at the sola r 
galactocentric distance of -Rq = 8 600pc l|Chen et al.ll200il ). 

The intrinsic I-band magnitudes of stars with spectral 
types from AO t o M9 were compiled f rom severa l sources 
in the literature llBessell fc Brettl l l988l : lGravlll992l : iLeggettI 
ll992l : lLegge"tt" et al .1 119961 ) while average colours were derived 
for the HST filters in a Section 3.1.2. We then populated a 
cone of length 7 200 pc in the direction of NGC 6611 accord- 
ing to the density behaviour described above, selecting stars 
randomly from the list previously generated. The opening 
angle of the cone simulates an area ten times larger than 



our real survey area to ensure that the contamination cor- 
rection plays no role in the statistical uncertainties of the 
IMF. 

The next step is to descri be the extinct i on pr ofile to- 
wards and behind NGC 6611. iBelikov et al] l|l999l ) shows 
that the reddening distribution of field stars towards the 
cluster has two peaks, one identified as a low-extinction fore- 
ground population and a second identified as a reddened 
background population obscured by the molecular cloud. 
The E{B — V) values measured towards both the massive 
cluster members and PMS stars seem to show a real spread 
in extinction. This suggests a first extinction screen associ- 
ated with the cluster's parent molecular cloud. The analysis 
in Section 3.1.1 strongly suggests that there is a second ex- 
tinction screen behind the cluster, probably associated with 
another spiral arm structure. 

We take all this information into account to construct a 
prescription for the behaviour of extinction as a function of 
distance along the line-of-sight towards NGC 6611. We as- 
sumed that E{B — V) has a double-step-like behaviour: it 
grows linearly with distance at the rate of 0.25 mag kpc~^, 
reaching Ay — 1.7 mag at the cluster position and then im- 
mediately increasing to Ay =3. 7 mag (in agreement with 
measurements towards cluster members). Behind the clus- 
ter, extinction continues to increase at 0.25 mag kpc~^ until 
the position of the second extinction screen. At a distance of 
2 500pc extinction climbs to Ay ~ 13mag, consistent wit h 
the 2MASS extinction profile from [Marshall et al. I (|2006l ). 
After that, extinction increases linearly to a maximum sim- 
ulated distance of 7 200pc. We adopt Ry = 3.75 mag up to 
and including the cluster position (the extinction behaviour 
is dominated by the average dust properties m the molec- 
ular cloud), while behind the cluster we adopt _Rv ~ 3.1, 
more consistent with the average properties of the Galactic 
interstellar medium (ISM). Each star in the simulation has 
extinction applied to its magnitudes according to its simu- 
lated distance and following this two-step function. Finally, 
the simulated apparent magnitudes are perturbed accord- 
ing to the observed distribution of measurement errors in 
the photometric catalogue. 

Fig. [7] shows a colour-magnitude diagram: filled and 
open circles represent the observed and simulated cata- 
logues respectively. The simulated catalogue seems to suc- 
cessfully reproduce the observed colour-colour and colour- 
magnitude diagrams in terms of the different stellar pop- 
ulations present: the foreground population, the moderate 
extinction background population and the heavily reddened 
population at larger distance behind the cluster. To test 
whether the simulated catalogue also quantitatively pro- 
duces the correct numbers of stars we investigated in more 
detail two particular regions in the IZ diagram: just below 
the PMS selection area and the region where the clump of 
stars with redder colours are observed. The comparison of 
the numbers of simulated and observed stars in colour bins 
is also shown in Fig.[7l For these two representative regions, 
the simulated catalogue produces stellar numbers that are 
entirely consistent with the observed catalogue. 

However, the simulated catalogue does not produce 
stars with more extreme red colours {I — Z > 2.8 mag). 
There are several possible reasons for this. The properties 
of the very reddened population are sensitive to the details 
of the extinction profile in the simulation, such as the lo- 
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Figure 7. Colour-magnitude diagram of the simulated and ob- 
served catalogues (open and filled circles respectively). Coloured 
symbols are the objects used for the comparison between the ob- 
served and simulated catalogues (respectively blue and red sym- 
bols). The bottom panel shows the I — Z colour histograms for 
two areas in the colour-magnitude diagram (dashed boxes): be- 
low the PMS selection region and at the location of the reddened 
star population. The simulated area is ten times larger than the 
HST field; the histogram has been renormalised to equal area. In 
terms of the number of objects in each colour bin the simulated 
(red histogram) and observed (blue histogram) catalogues agree 
very well. The simulation does not however produce objects red- 
der than about 1 — Z ~ 2.8mag (see text). Simulated objects that 
fall into the PMS selection region are potential contaminants, but 
are largely excluded on the basis of their reddening (see text). 

cation and magnitude of the second extinction jump, the 
fact that we do not account for any structure or dumpiness 
in the ISM and the adopted reddening coefficients. And at 
large extinctions, these coefflcients are manifestly not con- 
stant and they also depend on the intrinsic colours of each 
object. Furthermore, we only simulate main-sequence ob- 



jects; this reddened population has evolved stars associated 
with it, brighter and redder objects that we simply do not 
simulate. The properties of the first extinction jump affect 
only the population underneath the PMS selection region. 
Another source of uncertainty is the parameters eissociated 
wit h the Galactic ste llar density distribution, particularly 
/iR ICheneiaDliool). 

In spite of all these caveats, our simulation is able to re- 
produce the observed catalogue quite well. Fig. [7] shows that 
the number of possible contaminants (i.e. simulated objects 
that fall into the PMS selection region) is very small, even for 
such a deep survey. The heavily extincted background popu- 
lation is effectively removed by the J — H colour cut for PMS 
selection (Fig.[2}. By using reddening cuts (see Section 3.1.2) 
we are further able to exclude the majority of possible con- 
taminants from the foreground and moderately extincted 
background populations. The application of the reddening 
cuts also means that the exact details of the extinction be- 
hind the cluster are not crucial; as the area considered here 
is small, the number of foreground contaminants is so small 
that foreground extinction is also not an issue. However, 
some objects from the former population can have similar 
extinction as the young cluster, and thus contaminate the 
PMS selection region. These contaminants need to be taken 
into account when constructing the mass function. 

The procedure described here to mimic field star con- 
tamination produces a catalogue that is processed in exactly 
the same way as the observed stellar catalogue. A PMS and 
reddening selection is performed on the simulated catalogue: 
there are only 52 simulated contaminants in the PMS selec- 
tion region and of these only 17 or 14 fulfil the reddening 
requirements (respectively using the 2 or 3 Myr isochrones) . 
For these objects, intrinsic I-band magnitudes and stellar 
masses are estimated. The derived masses are used to com- 
pute a contamination mass function (MF). This MF is di- 
vided by 10 because the simulated area is ten times larger 
than our actual survey. We find that there are less than 2 
contaminants in the two lower mass bins. Even though con- 
tamination is low, we choose not to neglect it as it affects 
the more crucial lower-mass bins. The contamination MF is 
subtracted from the completeness-corrected IMF computed 
in the previous section (Section4.1). 

We consider other potential sources of contamination. 
Nearby giants are very bright and thus are saturated in the 
HST images. Giants associated with the heavily extincted 
population described above are not simulated. However, as it 
can be seen in Fig.[2J they are redder than J—H ~ 2 mag and 
thus do not contaminate the PMS selection region. Evolved 
stars between the cluster and this redder population could 
in principle contaminate the PMS selection; they are how- 
ever in negligible numbers due to the low stellar density and 
the small volume sampled. The reddening screen behind the 
cluster further precludes giants nearer the Galactic centre 
as well as bright galaxies from being PMS contaminants, by 
making these extremely red and faint. 

3.2 Higher and solar- mass stars 

We se lected from the IZJH catalogue of lOliveira et al.l 
(|2005l ) a subset of stars from the same area in the sky as 
the NICMOS mosaic in Fig.[T] No attempt was made to con- 
vert photometry from different filter systems to merge the 
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Figure 8. Colour-magnitude diagrams using tlie photometric catalogue from lOIiveira et alJ |200^. Small dots show the complete 
catalogue as described in that article and large circles represent the objects that fall onto the NICMOS mosaic (Fig. [T]l. The solid line s 
indicate the regions for selection of PMS candidates: filled circles show the PMS candidates. Solid and dashed lines are lSiess et al 
isochrones, respectively for 2 and 3 Myr. 



Table 2. NGC6611 higher-mass cluster candidates selected from the photometric catalogue of lOliveira et aL 1 1|2005|) . Positions are in the 
2MASS system. Columns provide the magnitudes and errors for /, Z, J, H, K and L' . E{B — V) = 0.70 mag is the extinction assumed 
for all PMS candidates. M2Myr ^iid MsMyr are the masses of PMS objects, assuming either an age of 2 or 3 Myr. The complete table is 
available in electronic form. 

ra dec / I — Z J H K L' A^2Myr -^^3Myr 

m sj (dm s-^ (j^^gj (j^^gj (j^^^gj (j^^^gj (j^^gj (j^^^gj (^^j (]y[^j 

18:18:40.09 -13:47:00.8 12.301 ±0.006 0.639 ±0.008 11.177 ±0.003 10.462 ±0.002 9.695 ±0.004 8.317 ±0.004 6.076 6.015 

18:18:40.61 -13:47:44.4 12.708 ±0.006 0.441 ±0.008 11.859 ±0.008 11.551 ±0.002 11.335 ±0.005 11.331 ±0.008 4.914 4.902 

18:18:42.25 -13:47:30.4 13.119 ±0.006 0.524 ±0.008 12.239 ±0.007 11.906 ±0.002 11.641 ±0.004 11.531 ±0.009 4.076 4.126 

18:18:39.37 -13:47:11.8 13.975 ±0.006 0.733 ±0.008 12.480 ±0.004 11.812 ±0.003 11.551 ±0.004 11.361 ±0.008 2.791 2.500 

18:18:39.83 -13:47:44.2 14.298 ±0.006 0.863 ±0.009 12.617 ±0.011 11.599 ±0.004 10.737 ±0.008 9.637 ± 0.007 2.707 2.420 

18:18:42.18 -13:47:22.4 14.845 ±0.007 0.868 ±0.009 13.042 ±0.006 12.193 ±0.002 11.832 ±0.005 11.593 ±0.009 2.474 2.285 



catalogues. This subset contains 176 objects from which we 
select PMS candidates based on colour-magnitude diagrams. 

3.2.1 Identification of PMS stars and membership 
selection 

As for the HST photometry, we follow criteria based on / 
magnitude, I — Z and I — J colours to select PMS objects 
from the groundbased photometry (Fig.|8} . J ~ H does not 
play a role in the PMS selection here as there is negligi- 
ble contamination from red stars in this magnitude range. 
The foreground field star sequence is scarcely populated in 
such a small area, indicating that contamination at brighter 
magnitudes is also very low. To increase the contrast of the 
field and PMS sequences, the diagrams also show the popula- 



tions for the larger cluster area as described in lOliveira et al.l 
(j 20Q5! ). Based on criteria described above, 150 stars are iden- 
tified as PMS candidates in the area covered by the NIC- 
MOS mosaic (Table[2l complete in electronic form only). 



3.2.2 Reddening determination of PMS stars 

No single set of isochrones covers the complete mass range 
covered by our HST and groun dbased catalogues. Fo r the 
magnitude range covered in the lOliveira et alJ l|2005h cat- 
alogue, the most appropr i ate se t of isochrones are those of 
ISiess. Dufour fc Forestinil l|2000l ) for solar metallicity that 
cover the mass range 0.1 — 7M0. This isochrone grid pro- 
vides m agnitudes calibrated accordin g to the conversion ta- 
ble from lKenvon fc Hartmamil l|l995h . Those magnitudes are 
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converted to the appropriate photometric system using pub- 
Hshed transformations: isochronal IR magnitude s were first 
conve rted to the 2MASS photomet ric system (ICarpenterl 
bOOll ) and then to the MKO system (|Tokunaga et al.ll2002l )~ 
I-band magnitudes need no conversion. 

We could have used the Siess isochrones and the IJH 
colour-magnitude diagram to determine the extinction to- 
wards the PMS candidates (theoretical Z magnitudes are 
not available), similar to what was described in Section. 3.1.2 
for the HST catalogue. However, and as expected for a 
brighter sample dominated by earlier type stars, the range of 
colours exhibited by these stars means that a large fraction 
of objects have a double reddening solution (see Sect. 3.1.2). 
In this regime, even small photometric errors can cause large 
uncertainties in E{B — V). Furthermore, we think there 
are unquantified systematic effects in the isochronal J and 
H magnitudes due to successive photometric transforma- 
tion. Taking all these factors into account, we consider this 
method to derive reddening rather unreliable for this par- 
ticular sample and therefore chose to adopt a constant red- 
dening E{B — V)= 0.7 mag (see Sect. 3.1.2). Using this ex- 
tinction value, intrinsic magnitudes are computed. Assum- 
ing an age of either 2 or SMyr and a distance of 1.8 kpc, 
we compute a stellar mass for each object using their intrin- 
sic I-band magnitudes and the lSiess et al.l l|2000l ) isochrones 
(Table 2)0 



3.2.3 Completeness correction and field star 
contamination 

Due to severe crowding towards the cluster centre and the 
relatively large number of saturated stars in the images, the 
groundbased data of this region are incomplete. To compute 
the completeness correction for the groundbased I-band cat- 
alogue we make use of the much deeper HST J and H cat- 
alogues. The I-band HST catalogue cannot be used in this 
way because it saturates at 7 ~ 16 mag, corresponding (very 
approximately) to 7 ~ 15.7 mag in the groundbased cata- 
logue. The procedure is rather simple. We start by assuming 
that the HST J and 77 catalogues are essentially complete. 
This is a reasonable assumption over most of the magni- 
tude range of interest here, as incompleteness in the IZJH 
HST catalogue is driven by the I- and Z-bands. The J-band 
HST catalogue saturates however at about J ~ 14 mag (very 
approximately J ~ 13.5 mag groundbased) thus we are un- 
able to study completeness at the very bright end. We use 
objects detected in both J and 77 and look for their po- 
sitional counterparts in the groundbased IZJH catalogue. 
This allows us to compute the recovery rate as a function of 
HST J-band magnitude. We need however an I-band com- 
pleteness function. Using an empirical linear approximation 
HST J-band magnitudes can be converted to MKO J-band 
magnitudes. And assuming a PMS age (Siess isochrone) the 
J-band magnitudes correspond uniquely to a groundbased 
I-band magnitude. The resulting I-band completeness func- 
tion for the groundbased catalogue has a similar shape to 
the HST I-band completeness function shown in Fig. [S] but 
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6 lOIiveira et al] | |2005|) adopted a distance of 2 kpc, Ay = 3.4 mag 
{E{B — V) = 0.9 mag) and an age of 3Myr. These values are 
different from those adopted here so estimated masses will differ. 
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Figure 9. Completeness and contamination corrections to be ap- 
plied to the computed IMFs, for the HST (solid line) and ground- 
based (dashed line) catalogues. The histograms show mass func- 
tions (in black), corrected for incompleteness (red) and contam- 
ination (blue) for the HST catalogue only (see text). Contami- 
nation is small throughout, even in the brown dwarf mass bins. 
These histograms refer to the 2 Myr calculations, similar ones 
were constructed using 3 Myr isochrones. 



shifted to brighter magnitudes: it has a value just below 
unity to about 7 ~ 15 mag, falling after that as a second 
degree polynomial to a value of about 0.1 at 7 ~ 19 mag. 
To construct the IMF, each PMS candidate is ascribed a 
weight equal to the inverse of the completeness correction 
at its observed I-band magnitude. 

Field star contaminatio n in the ma gnitude range of in- 
terest here is very small. iQliveira et al.l (2005) used IZ pho- 
tometry of an offset field to estimate foreground contam- 
ination. Following the same procedure, we estimate that 
~ 2 objects might be field star contaminants in the NIC- 
MOS mosaic area in this magnitude range; this is an upper 
limit as we use only IZ criteria (no offset J77 photometry 
is available). Such low level of contamination is consistent 
with what we found for the HST catalogue (Section 3.1.4). 
We chose not to apply any contamination correction to this 
data set, as we have no constraints on the likely masses of 
possible contaminants. 



4 THE INITIAL MASS FUNCTION OF 
NGC6611 

4.1 Constructing the cluster IMF 

In the previous section, we identified PMS candidates in 
NGC6611 and computed th eir intrin sic I-band mag nitudes. 
Using isochrones from Bar affe et al. I 11998) and Siess et al.l 
(|200C[ ) (respectively for the HST and groundbased datasets) 
we computed masses assuming a distance of 1.8 kpc and us- 
ing either 2 or 3 Myr isochrones (Tables[T] and [2} . We also 
estimate the contributions from completeness and contami- 
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Figure 10. IMFs derived assuming an isochronal ages of 2 Myr (left) and 3 Myr (right). IMFs were constructed independently using 
the HST (blue) and groundbased (red) catalogues. The extreme bins in each histogram should be disregarded because they are not 
completely sampled by the catalogues. The two contributions agree remarkably well. The choice of isochrones also does not seem to 
influence the derived IMF in any significant way. 



nation. Fig. [5] illustrates the corrections that are applied to 
the observed IMF to account for these two effects (in this 
particular case for the 2 Myr set of calculations). Firstly, his- 
tograms are constructed for the PMS masses (bins of 0.3 dex 
in log mass) : with unity weights (black histograms) and with 
weights equal to the inverse of the completeness function 
(red histograms) to perform the completeness correction. As 
described in Sect. 3.1.4, a contamination mass histogram is 
computed and is then subtracted from the completeness- 
corrected histogram (blue line). The contamination level is 
low throughout; in the two lower mass bins there are re- 
spectively ~0.2 and ~ 1.2—1.5 contaminants (less than 10% 
contamination). The final histograms (i.e. with complete- 
ness and contamination corrections applied) are the IMFs 
to be discussed in the following paragraphs. We apply no 
correction to account for unresolved binaries. The bin size 
of 0.3 dex was chosen to match other published IMFs (Sec- 
tion 4.3); however we tested that this choice does not sig- 
nificantly change the measured properties of the IMF. The 
extreme bins in each histogram should be ignored because 
the catalogues do not fully sample the mass range in those 
bins. 

The groundbased and HST samples overlap in 
the mass range 0.15 1.5 Mq. The resulting IMFs 
are shown in Fig. 1101 (log N versus log M /Mq ) , as- 
suming an age of 2 or 3 Myr. There is a very good 
agreement between the derived IMFs, for the bins 
where they overlap, despite different datasets and 
photometric systems, reddening treatments, theo- 
retical models etc. For these two mass bins we aver- 
age the two contributions to obtain the final IMFs. 

Fig.[lO]also shows that the choice of isochronal age does 
not make a significant difference in the constructed IMF 
properties. The most noticeable difference is the somewhat 
sharper fall towards higher masses for the 3 Myr determina- 



tion. However, the two histograms are consistent within the 
statistical errors and furthermore these higher mass bins are 
affected by low-number statistics and sampling effects. As 
the two IMFs are indistinguishable, we settle on the 2 Myr 
representation for further analysis and only when relevant 
provide comments using the 3 Myr-derived IMF. 

We find that the IMF in NGC 6611 rises Salpeter-like at 
higher masses, consistent w ith many determina tions avail- 
able in the literature (e.g., iDufton et all 12006*). The IMF 
flattens out between ~ 1.25 and 0.3 Mq (with the most likely 
peak position between 0.3 and 0.6 Mq) and then drops into 
the brown dwarf regime. There is an apparent secondary 
peak at ~ 0.05 Mq. However, as we disregard the lower mass 
bin, there is no evidence for a steep decline in the IMF to- 
wards lower masses. 

We identified 30—35 brown dwarfs (respectively assum- 
ing 3 or 2 Myr) in our observed sample (before any complete- 
ness correction was applied, see Table[l|. The brown dwarf 
to star ratio, defined as the ratio of the number of stars 
with masses iii the r anges 0.02-0.08 Mq and 0.08-10. Mq 
(|Bricefio et al.ll2002l ). is 0.19 ± 0.03 (0.17 ± 0.03 for 3 Myr) 
The star to brown dwarf ratio, defined by lAndersen et al.l 
(2008) as the ratio of the number of stars with masses in the 
ranges 0.08-1. Mq and 0.03-0.08 Mq is 4.4i;;;^ {^-itol for 
3 Myr). In the next subsections we will compare these ob- 
served properties both with theoretical IMF representations 
and with other young clusters. 



4.2 Comparison with IMF theoretical 
expectations 

The most commonly used theoretical representations of the 
field IMF ma ke use of a multiple-part power-law (e.g., 
Kroupal [200 ll) or a Scalo-like log-normal function (e.g., 
Chabriet 20051 ). Both functional forms describe the IMF as 
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Figure 11. Comparison of the IMF of NGC 6611 witli well known 
theoretical funct ional forms: th e multiple-part power-law IMF (la- 
belled bl ue line. |Kroupal[200lh and the log-normal IMF (labelled 
red line, IChabrieJf2005r ) . We have fitted such functions to the 
observed IMF (ignoring the two extreme bins); the most ade- 
quate fits are overplotted on the IMF. Both when compared to 
the power-law and log-normal formulations, the NGC 6611 IMF 
peaks at a slightly higher mass, ~O.45M0 (see text). 



reaching a peak at a few tenths of a solar mass and falling 
both into the higher-mass and brown dwarf regimes. The 
iKroupal (|200ll ) power-law IMF (in linear units) is charac- 
terised by the index a and is defined as ^{M) = dN/dM oc 
M~"' . The average Galactic-field IMF (uncorrected for un- 
resolved binaries) is Salpeter-like to 0.5 Mq with qi ~ 4-2.3 
— in this representation the Salpeter slope is equal to -1-2.35; 
between 0.08 and 0.5 Mq, the IMF flattens somewhat, con- 
tinuing to rise but less steeply with 02 ~ +1.3; t he IMF 
drop s into the brown dwarf regime with Q3 ~ +0.3 (jKroupal 
I2OOII ). The log- normal description (in logarithmic units) has 
the form f = d N/dlogM oc exp(-(log m - log mc)^/2cr^). 
IChabried (|2005') find that mc = 0.25MQ and a = 0.55 are 
required to reproduce the stellar MF in the solar neigh- 
bourhood (no correction for unresolved binaries). This log- 
normal distribution is also found to fi t well the MFs of a 
number of open c lusters like Blanco 1 (|Moraux et al.l [20071 ) 
and the Pleiades (|Moraux et al.ll2003l ). 

In Fig.[Tl]we compare the IMF of NGC 6611 with both 
the power-law and log-normal IMFs. We stress that we 
have applied no correction to account for unresolved binary 
systems. In terms of a power-law representation, we find 
that the best fits are provided by power-law indices a; of 
+2.36, +1.18 and +0.52, respectively for masses > 0.9 Mq, 
between ~0.9 and ~0.45Mq and < 0.45 Mf^^. These iii - 
dices are consistent with those determined bv lKroupal (|200ll ) 
within the uncertainties. However, the IMF in NGC 6611 
flattens at higher masses (~0.9Mq instead of 0.5 Mq) and 
starts to fall sharply also at higher masses (~ 0.45 Mq in- 
stead of 0.1 Mq). When using the log-normal function we 
determine rric - 0.40 ± 0.04 Mq and a = 0.56 ± 0.04. While 



the width of the distributions is similar in NGC 6611 and 
in the field, the peak occurs at higher masses in NGC 6611. 
Both the comparison with the power-law and log-normal 
IMFs suggests that the IMF in NGC 6611 peaks at higher 



4.3 Comparison with other young clusters 

In this section we compare the IMF of NGC 6611 with those 
of other young clu sters and assoc iations com piled from the 
literature; Taurus j|Luhman|[2003 ). th e ONC |Muench et aD 
|2002| ; ISlesnick et all |2004^ IC348 (iLuhman et al.l |2003|) 
ChamaeleonI llLuhmarill2007l ) and NGC 2024 (jLevine et all 
liool). ISlesnick et all (l2004l) construct e d the IMF for the 
ONC using the iD'Antona fc MazzitellH (| 19971 ) evolutionary 
tracks ; for the remaining determinations the iBaraffe et al.l 
(Il998l ) models were used. These clusters are younger than 
a few Myrs in order to minimise efi^ects like dynamical evo- 
lution and mass segregation. No correction for binarity 
was applied to the IMFs of any of the clusters dis- 
cussed here, assuring that these IMFs are compara- 
ble, so long as their binary properties are not too 
differentO 

In Fig.[l2]we show the IMFs of the clusters considered 
here. A noticeable difference is that the IMFs of both Taurus 
and NGC 6611 (Fig.[l2l left) seem to peak at slightly higher 
masses than those of the other young clusters (Fig. 1121 right); 
the IMF peaks are at ~0.8Mq and ~0.5Mq, respectively 
for Taurus and NGC 6611 and at O.l - 0.2 Mq for the other 
clusters. Another quantity also used to characterise the IMF 

— and that might be less dependent on the adopted binning 

— is the characteristic mass, defined as the mid-point of the 
IMF plateau (Elmegreen cJ^aL 2008); for most clusters in 
Fig.[l2]it occurs at ~0.3Mq while for NGC 6611 it is at 
~0.45 Mq. 

We perform a Kolmogorov-Smirnov statistical compar- 
ison to test whether the observed IMF differences are signif- 
icant. When the IMF of NGC 6611 is compared to those of 
Taurus and the ONC, the probability that they are drawn 
for the same cumulative distribution is 0.1 and 4 x 10~^, 
respectively. The probability that the IMFs of the ONC and 
IC 348 are drawn from the same distribution is 0.65, confirm- 
ing that indeed they are very similar. This analysis confirms 
that the shape of the IMF of NGC 6611 more resembles that 
of Taurus and it is very unlikely to mirror that of the ONC. 

We can also compare the relative numbers of stars and 
brown dwarfs for these clusters. We should point out that 
such ratios are problematic to compute as it is difficult to 
achieve survey completeness; for instance Tau r us is a very 
spatially extended region and ICambresv et al.l (|2006h claim 
that current IC 348 surveys are still incomplete. We list 
available ratios for the 6 clusters in Tabled The brown dwarf 
to star ratio varies from 0.13 for IC 348 to 0.30 for the ONC 
and NGC 2024. The brown dwarf to star ratios of the ONC, 
NGC 2024 and Chamaeleon I are similar within the uncer- 
tainties while those for Taurus and IC 348 are lower. The 



S IChabrierl (l2005t) has shown that when unresolved bi- 
naries are taken into account the peak of the disc 
IMF moves to sHghtly lower masses (from 0.25 Mq to 
0.20 Mq). 
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Figure 12. Comp arison between IMF s for several young clusters and associations: NGC 661 1 (this work, red), T aurus l lLuhmarJl2004l 
green), the ONC llMuench et al.ll200l orange soli d line, rebinned to match other cluste rs: ISlesnick et al lj2004l . orange dashed line), 
IC348 llLuhman et al.ll2003l. "cyan). ChamaeleonI llLuhmanI |2007| . blue) and NGC 2024 llLevine et al.ll2006l . maeenta). The y-axis is 
in arbitrary log A'^ units and published cluster IMFs have been renormalised for display purposes. The vertical dotted lines indicate the 
typical positions of the plateau and the characteristic mass Mc for most clusters. The peaks of the IMFs of Taurus and NGC 6611 seem 
to occur at higher masses. 



ratio we determine for NGC 6611 is consistent with that of 
Taurus. The ratio of stars to brown dwarfs (as defined by 
lAndersen et al.l l2008h varies from S.StES:? for the ONC to 
8.3l^:^ for IC348. The ratio for NGC 6611 is between those 
for ChamaeleonI and Taurus. 

An issue still being debated is whether the shape of the 
IMF shows any imprint of the conditions under which the 
star formation process takes place. Those initial conditions 
are impossible to constrain directly, but we can try to infer 
them from the observed present-day properties of the result- 
ing stellar populations. We use cluster data available in the 
literature (|Andersen et al.l l2008'. and references therein) to 
crudely estimate an average stellar surface density of PMS 
stars (dividing the number of stars with masses in the range 
0.03— IMq by the area) for each cluster. We also estimate 
the average ionising flux at the position of the PMS stars 
for NGC 6611 and the ONC — using tabu lated total ionising 
fluxes from ,Osterbrock fc Ferlandl (|2006| ): the other popula- 
tions contain no O-stars. These and other relevant cluster 
properties are listed in Tabled 

The ONC is the most extreme of these associations 
in terms of present day stellar density. ChamaeleonI, 
NGC 2024 and IC 348 all have similar stellar densities and 
are less dense than NGC 6611 and the ONC. Taurus on the 
other hand is the prototypical example of a low density 
quiescent star forming environment. The average ionising 
fluxes for NGC 6611 and the ONC are similar, despite the 
fact that NGC 6611 has more and more powerful ionising 
sources. NG C 6611 is the most m assive of the clusters con- 
sidered here (jBonatto et al.ll2006l ). 

The IMFs of the ONC, ChamaeleonI and IC348 are 
similar, despite their different stellar densities; although 
compared to these other clusters IC 348 has an abnormally 



low number of brown dwarfs. The IMF of NGC 6611 most 
resembles that of Taurus, both in terms of peak mass and 
relative number of brown dwarfs, but the star formation en- 
vironments of NGC 6611 and Taurus are dramatically differ- 
ent. The IMFs of NGC 6611 and the ONC are very different, 
despite both being examples of low-mass star formation ac- 
companied by O- and B-type stars. 

The properties described in the Tableware present day 
cluster properties and the ages of the clusters span 1 — 3 Myr. 
Some dynamical evolution has probably occurred so the pop- 
ulations were likely denser and, in the cases of NGC 6611 and 
the ONC, the PMS stars were likely nearer to the ionising 
sources. We only analyse a small central region of NGC 6611 
and it is not possible to assess at the moment whether the 
IMF we constructed is representative of the whole cluster. 
Still, no clear trend has emerged so far that relates IMF ob- 
served properties with known cluster properties like ionising 
radiation or stellar density. Our results suggest that the ra- 
diation field from the O-stars does not play a dominant role 
in shaping the lower-mass IMF. 

4.4 Discussion and implications 

The goal of this paper was to investigate the properties of 
the low-mass IMF of NGC 6611, and in particular the brown 
dwarf frequency and peak or characteristic mass. The intrin- 
sic complexity of the IMF and the still limited understanding 
of which processes are dominant in shaping it (iBonnell et al.l 
2007; Elmcgreen 20Q§) means that it is not possible to model 
an individual IMF (except by using analytic parameterisa- 
tions like the ones described in Section 4.2). Instead by com- 
paring the IMFs of clusters and associations with different 
present-day properties we can try to constrain which phys- 
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Table 3. Compiled cluster and IMF properties (relevant t o th e c omparison descr i bed in S ection 4.3). R e ferenc es for the ratios (see 
text for definit i ons) a re as follows : 1 thi s w ork; 2[L uhman ( 2000l) : 3 lAndersen et al.l ([Sooi); 4 lLevine et~al] | |2006|): siL uhman (2007); 6 
iLuhman et al.l l|2003l ): TjLuhman' (2004); 8 Guiou ot al. (20c 3). The other cluster properties ar e compiled from Andersen ct al. (200p, 
and references therein) ; iLuhman (2008) ; Hcrbst. (2008) ; Meyer et al.l l l2008t) ; iMuench et al] | |2008|) ; lOliveiral |2008) ; Kcnyon ct al. (20oi). 
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ical conditions and eventually which physical process are 
more important. 

As recently reviewed by lElmegreenI (|2008l ). there are 
three main types of star formation theory: fragmentation 
(both gravitational and turbulent), competitive accretion 
and interruption (that includes both ejection and photo- 
evaporation scenarios). These theories need to be able to 
explain not only the properties of the IMF, but the binary 
and circumstellar disc fractions across all masses and the 
dynamical properties of stellar populations. 

There is abundant observational evidence that sup- 
ports a similar forma tion mechanism for brown dwarfs 
and solar- mass stars ()Luhman et al.l l2006h . However, an 
open question in the understanding of brown dwarf (and 
low-mass star) formation is how the low-mass cores are 
stopped from accreting further so that their final mass can 
remain sub-stellar. Two mechanisms have been proposed 
to halt accretion: ejection of newly formed fragments and 
phot o-evaporation of collapsing cores. The ejection scenario 
(e.g., iReipurth fc Clarke! I2OO1I ) proposes that the smallest 
clumps in multiple systems are ejected and are thus de- 
nied access to the molecular gas reservoir. There are sev- 
eral problems with this scenario, namely that the spa- 
tial distribution of solar-ma ss stars and b r own dwarfs in 
clust ers seems to be similar (jLuhman 2004'; Luh man et al.l 
I2OO6I ). and that both the binary (iBurgasser et al.ii2007l ') and 
^properties J .Tavawardhana et al. l l200,- j) of the two pop- 
ulations are a lso similar. The photo-evaporation scenario 
l|Whitworth fc Zinnccke] |2004| ) relies on the proximity to 
massive O-stars to stop accretion. This would imply that 
we should find a higher brown dwarf to star ratio in OB as- 
sociations than in quieter regions; this is not supported by 
our observations in NGC 6611. 

Recently iBonnell et al.l ()2008h suggested that gravita- 
tional fragmentation of infalling gas into stellar clusters can 
form low-mass stars and brown dwarfs without invoking any 
additional mechanism. Increased gas density within the col- 
lapse region gives rise to a filamentary-like structure within 
which the Jeans mass is lower, allowing lower-mass clumps 
to form; the infall velocity of the gas prevents the low-mass 
clump from accreting significant amounts of gas, preserving 



its low-mass or brown dwarf status. Another mechanism like 
turbulent fragmentation would howeve r be necessary to form 
brown dwarfs in low-density regions l|Bonnell et al.|[2008l ). 
As discussed in the previous section, the observations do not 
suggest any variation of brown dwarf fraction with present 
day average stellar density. 

The fact that massive star feedback does not seem 
to have a significant impact on low mass star forma- 
tion is not necessarily surprising. The idea that compet- 
itive accretion in a cluster environment is an important 
process in the formation of the m o re massive stars is 
gaining favour (iBonnell fc Bate! l2006l : iBonnell et al] |2007| ; 
lElmegreenll2008l ). In this context, fragmentation produces 
lower-mass stars and a few higher-mass stars form subse- 
quently by accreting more mass in the cluster gravitational 
potential. This would imply that by the time the most mas- 
sive stars are fo rmed, the main properties of the IMF ar e 
already defined (|Bonnell fc Batell2006l ; IBonnell et al.ll2007l ); 
we would then not expect significant differences between the 
IMFs of OB associations and quiet star forming regions, con- 
sistent with what was shown in Section 4.3. 

The shape of the IMF between 0.1 and IMq also 
presents a challenge. The characteristic mass of the IMF 
Mc is observed to be essen tially constant for most star form- 
ing regions (see review by Hi megreen et al ] l2008l ). Based on 
numerical simulations it has been proposed that this char- 
acteristic mass (and the associated plateau) is related to 
the therm al Jeans mass at the onset of (isothermal) col- 
lapse (e.g.. lBate fc Bonnellll2005l ). One could then intuitively 
expect the plateau properties to vary with the environ- 
mental conditions, via for instance the dependence of the 
Jeans mass on the den sity in the molecular core. However, 
lElmegreen et al] (|2008l ) show that, if grain-gas coupling is 
taken into account, the thermal Jeans mass depends only 
weakly on environmental factors like density, temperature, 
metallicity and ra diation field. Indeed, recent simulations by 
Bon nell. Clarke fc Ba te (2006), using a more realistic equa- 
tion of state to allow for the coupling of the gas and dust 
at the high densities typical of star forming regions, predict 
a characteristic mass that is relatively independent of the 
initial conditions for star formation. 
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Fig. [12] shows some significant differences in tlie ob- 
served IMFs of young clusters. At the present, and in the 
absence of any trend of IMF properties with cluster proper- 
ties, we remain unable to explain such differences. It is pos- 
sible that the lower-mass IMF is shaped by a combination of 
different processes and that a particular effect might become 
dominant depending on the environmental conditions. 



5 SUMMARY 

We present results based on the deepest images to date in the 
central region of NGC 6611. The observations were obtained 
with ACS/WFC and NICMOS on board HST and reach 
down to / ~ 26 mag. We use this photometric catalogue 
to construct the cluster IMF from 1.5 MQ_to_a02 Mq. We 
use the photometric catalogue from iOliveira et al.l |200^) to 
extend this IMF to higher masses (~7M0). 

We find that the IMF of NGC 6611 is well described by 
a lo g-norma l distr ibution. However, both w hen comp a red to 
the lKroupal l|200ll ') multiple power- law and lChabrieil (|2005| ) 
log-normal distributions, the IMF of NGC 6611 seems to 
peak at higher masses than the Galactic field and the solar 
neighbourhood IMF. The higher-mass slope is Salpeter-like 
and we find no evidence of a sharp decline towards the brown 
dwarf regime. We identified 30 — 35 brown dwarf candidates, 
depending on the assumed cluster age. 

We compare the IMF we derive for NGC 6611 to those of 
other clusters available in the literature. We find conclusive 
evidence that the IMF of NGC 6611 more closely resembles 
that of Taurus than that of the ONC. This is surprising 
since Taurus is the prototypical low-density star-forming re- 
gion while both NGC 6611 and the ONC are examples of 
low-mass star formation in OB associations. Even though 
the present day properties of these two clusters (stellar den- 
sity and ionising radiation field) are similar, a K-S test shows 
that their observed IMFs are not drawn from the same mass 
distribution. Our analysis yields no trend that relates either 
present-day stellar density or intensity of the ionisation field 
with observed IMF properties. Still, the fact that the IMFs 
of Taurus and NGC 6611 are similar suggests that the forma- 
tion of lower-mass stars can be unaffected by their massive 
siblings. 

Even though there is mounting evidence that suggests 
variations in the lower-mass IMF for young clusters, we are 
still at present unable to fully explain them. A large sam- 
ple of clusters, observed and analysed consistently, is neces- 
sary to try to investigate environmental links in any detail. 
New instrumentation is already allowing us to investigate 
the lower-mass IMFs in extragalactic environments. In that 
context IMFs constructed solely on the basis of photometric 
catalogues will become increasingly important, as they offer 
the only reliable comparisons with extragalactic determina- 
tions. 
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